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1. Introduction. The most common method for estimation of the lower eigen- 
values of a differential operator is the Rayleigh-Ritz method, or a variant such as finite 
elements in engineering problems, or Hartree-Fock in quantum mechanical problems. 
The emphasis in many numerical studies is on obtaining accurate eigenvalue approxi- 
mations in an efhcient and cost effective manner. Herein, we explore a different ques- 
tion, namely, is the calculation useful if needed to estimate a bifurcation coefficient in 
a nonlinear branching problem? If the eigenfunctions of the corresponding linearized 
problem must be approximated, the strength of the nonlinearity must be considered. 
In particular, the eigenfunctions of the linearized problem need to be approximated 
in a topology strong enough that nonlinear quantities to be calculated behave contin- 
uously. This may or may not be the case with the standard Rayleigh-Ritz method. 
We know of one prior paper on this question [S] , in which the nonlinearity is taken to 
be continuous in the underlying Hilbert space topology. 

Our approach is to introduce a whole scale of Rayleigh quotients using real pow- 
ers of the selfadjoint operator which determines the corresponding linearized problem. 
This is a version of quadratic interpolation, also called Banach space interpolation be- 
tween Hilbert spaces. A discrete set of such Rayleigh quotients was called "Schwartz 
quotients" in [16] , and used to generate a theoretical algorithm for eigenvalue estima- 
tion that is closely related to the power method for matrices. We use our somewhat 
esoteric construction to seek sharp convergence rate estimates for bifurcation coeffi- 
cients. This proceeds by use of eigenfunction estimates in Hilbert space topologies 
strong enough for the nonlinear approximation procedure to converge. 

In the next section we describe the basic Rayleigh-Ritz method, and formulate a 
model problem concerning a nonlinear rotating string with variable density. In Sec- 
tion 3 the scale of Rayleigh quotients is presented, and relevant facts about quadratic 
interpolation between closed subspaces of Sobolev spaces determined by homogeneous 
boundary conditions are stated. The basic eigenvalue-eigenvector convergence theo- 
rem is also presented here, with the proof attached as the appendix. Then in Section 
4 we prove theorems on spectral implementation of the basic convergence theorem. 
This method is known as the spectral Galerkin method in the numerical analysis 
literature (cf. [10], [14], [28]). Most of this section can be implemented via finite 
elements, which we comment on in Section 6. In Section 5 we analyze the nonlinear 
rotating string problem, and present a few numerical results for this problem with a 
particular density function. Finally in Section 6 we illustrate our theory in the context 
of bifurcation problems for some partial differential equations, and close with a few 
remarks. 

2. The Rayleigh-Ritz Method and a Model NonUnear Problem. Let H 

be a separable complex Hilbert space with norm || • || and inner product (•,■). Let A be 
a selfadjoint operator in H with domain V^A). We assume that A is positive definite, 
i.e., \\Au\\ > a||u||, u S 'D{A), a > 0, and that the lower portion of the spectrum of 
A consists of isolated eigenvalues of finite multiplicity. Thus the lower spectrum of A 
may be written 

< a < Ai < A2 < . . . < Aoo 

with corresponding orthonormal eigenvectors iti, U2, ■ • ■ , where Aoo denotes the least 
point of the essential spectrum of A. By convention, Aqo = cxd if A has compact 
inverse. Further, let a{u) be the quadratic form corresponding to A, i.e., the closure 
of {Au,u), with domain I?(a). Then a{u) > for u G 'D{a), and we let a{u,v), 

u,v ^ 'D{a), be the corresponding Hermitian symmetric bilinear form. 
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The Rayleigh-Ritz method for estimation of eigenpairs of A begins with a family 
of trial vectors {wi}2^i C I?(a). It then proceeds with the setup and resolution of the 
generalized matrix eigenvalue problem 

[a{wt,Wj)] X = A[{wi,Wj)] X. 

If the resulting matrix eigenvalues are denoted by A^"^ < A^"^ < . . . < a1"\ the first 
monotonicity principle 37 guarantees that for each i — 1,2, . . . ,n, Xi < A[^\ thus 
providing upper bounds for the lowest n eigenvalues of A. Let u']^\u'"^\ . . . be 
the corresponding orthogonalized eigenvectors, normalized in Ti.. Thus {u[^\u'j^^) = 

5ij and a(u["\ ^^j"^) = ^1"'' where Sij is the Kronecker symbol. 

Completeness of {wij^^i in I?(a) is sufficient to guarantee that for each eigenvalue 
Xi of A below Aoo, a|"^ Xi and — > Ui in X'(a) both as n ^ cxd. (cf. [4] [15] 
[37]). Rate of convergence estimates are also available from these sources. But to 
use the Rayleigh-Ritz eigenvectors for other calculations of interest requires estimates 
in other norms, which is the thrust of this paper. Herein we look specifically at the 
estimation of bifurcation coefficients in nonlinear eigenvalue problems. Convergence 
in the energy norm, i.e., that of I?(a), may be insufficient to handle the nonlinearity, 
or when sufficient may give only a crude estimate. To illustrate this point we now 
begin the examination of a nonlinear rotating string problem. 

Perhaps the simplest nonlinear refinement of the standard linear model for a 
tightly stretched flexible string of length L with linear density p, rotating with uniform 
angular velocity uj about its equilibrium position along the x— axis is 

Herein, y is the deflection from equilibrium at the point x, and the condition of force 
equilibrium in the x direction (in the absence of external forces in that direction) 
implies that the tension T satisfies 

dx 

T —— = Tq = constant, 
as 

so that by solving for the deflection y, the tension T is given by 

(cf. [23]). As usual, s denotes arclength. When p is constant, (|2.1[) is solvable 
explicitly in terms of elliptic integrals [19j . In this section we examine calculation of 
small, but not inflnitesimal, deflections when p is variable. Setting A = lu'^/Tq and 
y — \ftu, with e small and positive, yields the nonlinear eigenvalue problem 



— u" = Xp(x)u \/l + eu''^ ~ 0, ' — 

dx 

and we will require flxed end conditions 



u(a) — u(b) = 0, 
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signifying that the equihbrium position of the string is the interval [a, b] of length L. 
We assume that p is smooth and strictly positive on [a, b]. 

To better illustrate applicability of the theory we will in fact consider the family 
of nonlinear eigenvalue problems 

- m" = A p{x)u (l + e {u'f^ " = 0, u{a) ^ u{b) ^ 0, (2.2) 

where i is a positive integer. These are all W^'^{a, b) perturbations of the linearized 
problem obtained by setting e = 0, and the rotating string model is given when t = 1. 

Solutions of (|2.2p are smooth and the eigenvalues of the linearized problem are 
simple, so known results of bifurcation theory employing the Banach space C^[a, 6] 
imply that there are solution pairs (AcUe) analytic in e for small e emanating from 
any eigenpair (Ao,mo) of 

-iL" = \p{x)u, M(a) = u(6)=0, (2.3) 
(cf. [3T], [Mj)- So we write solutions of (|2.2p as 



(2.4) 



where — u\^=Q^ vi — ^|e=o, etc. As we shall see, calculation of the coefficients 
vi,iyi, etc. in (|2.4p based on Rayleigh-Ritz type solutions of (|2.3p will require Ritz 
methods convergent in norms stronger than that of the energy norm, i.e. the norm of 
i/p (a, b). That this is the case is readily seen by expanding the nonlinearity in (|2.2p 
to obtain 

// N / 1 , ^ / /\2t , (1 ~ 2t) 2 / M4t , 



Successive equations for j/q, wq, i^i, t^i, ■ ■ • are obtained, as usual, by equating coeffi- 
cients in e in this expression. The first equation obtained is the linear eigenvalue 
problem 

- = p(x)vt), va{a) = vo{b) = 0. (2.5) 

Our goal is to investigate the use of Rayleigh-Ritz methods for (|2.5p in the calculation 
of the bifurcation coefficient vi , and further terms in the two series (j2.4p . The second 
equation is 

- v'l - 1^0 p{x)vi = 1^1 p{x) va + ^ p{x) vo {vQf\ vi{a) =^ vi{b) ^ 0. (2.6) 

The Fredholm alternative yields solvability of (|2.6[) if and only if 

1^1 {p{x)vo,vo) + ^ {p{x)vo{v'o)'^\vo) = 0, 

where (• , •) is the usual L'^{a, b) inner product with respect to Lebesgue measure dx. 
Taking vq to be normalized in with respect to the measure p{x)dx over (a, b) (which 
we will take as the space Ti with inner product (•, •)), i.e., 

{p{x)vQ,vo) ^ {vq.vq) ^ I, (2.7) 
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gives z/i as 

^1 = -| {p{x)v,{vX,v,) = (^;oK)'',^;o). (2.8) 

The component of vi in the "p orthogonal complement oivo^ is uniquely determined 
by solving (|2.6p with this value of I'l. Note that I'l is negative, i.e., the bifurcation is 
subcritical. By also setting the component of vi parallel to wq in Ti. equal to zero, one 
obtains that the parameter e^t is the Ti norm of the component of the power series 
solution of (|2.ip in the direction of wq. 

The third equation obtained from (|2.4|) is 



-W2 - Vq p{x)v2 = V2 p{x)vo + 2uip{x)vi + y p{x) Vq (v'^f* 

+ y P{x) vi {v'of* + 2vo p[x) vq v[ (wq)^*"^ 

i^o{l -2t) I it ,„ „s 

+ p{x)vo{vq) , (2.9) 

with V2{a) = V2{b) = 0. Retaining the normalization (|2.7|) . the Fredholm alternative 
yields solvability of (|2.9p if and only if 

V2 = -2l/i(ui,Wo) - y (wo (Wo)^* ,Wo) 

- y (^^1 ("o)^* ' ""o) - 2i^Q {vo v[ {v'^f*'^ , vo) 

i^o (Vo) ,vo). (2.10) 

With this value of 1^2, V2 is uniquely determined in the Ti. orthogonal complement of 
vq, while in the direction of ?;o the component of V2 is taken to be zero. This preserves 
the previous interpretation of the parameter £2* . 

Note that if (|2.5p is solved approximately by the Rayleigh-Ritz method, conver- 
gence of the approximate eigenvectors is in the topology of i?o (a, b). This is sufficient 
to guarantee convergence of the corresponding approximation to i>i as calculated from 
(|2.8p only for t = 1. To rigorously approximate z/i for t > 1, or 1^2 for any value of 
t, requires that approximate eigenvectors converge in CQ[a,6], which is strictly con- 
tained in HQ{a,b). We will obtain such convergence estimates by using quadratic 
interpolation between Hilbert spaces to generalize the venerable notion of Rayleigh 
quotient. This will yield convergence rates in H°'{a,b) with a > 3/2, which implies 
convergence in Cq [a, b] via basic Sobolev theory. Then rates of convergence of these 
nonlinear functionals will be estimated. 

3. Quadratic Interpolation and Fractional Rayleigh-Ritz. To provide a 
theoretical framework for the eigenvector estimates required in the preceding example, 
let T = A~^, which is a bounded selfadjoint operator in Ti. It is well known that 
T is also selfadjoint in I?(a) with the topology induced by the inner product a{u,v) 
i.e., the restriction of T to 'D{a) is bounded and satisfies a{Tu, v) = a{u, Tv) = (u, w), 
li, w, G 'D{a). The Rayleigh quotient for reciprocals of eigenvalues of A is 

{u,u) a(Tu,u) 
a{u) a(u) 

In this form the Rayleigh-Ritz method produces lower bounds, 

-1 
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as n — > oo from trial vectors {wi}°^^ that are complete in I?(a). Now for real r 
let be the r**^ power of A as defined by use of the spectral theorem [21]. This 
is particularly simple to describe if A has a complete set of eigenvectors {wij^i, 
orthonormal in H. Then u = otiUi G Ti. if and only if J^ili l^^iP < ^^'^ if 

T > 0, It e V{A^) if and only if J^Zi < while for r < we take X>(yl^) 

to be the completion of H, in the norm jj^'^ulp = X^i^i |c«ip. Herein a; = (u, u;). 
Whether or not A has a complete set of eigenvectors, 'D{A^) and 'D{A~'^) are dual 
spaces, where duality is taken relative to (•,•). Since A has a strictly positive lower 
bound, if A is bounded all of these norms generate the same topology. Our interest is 
in unbounded A, in particular differential operators, and then if ti > T2, the topology 
corresponding to ti is strictly stronger than that for T2 ■ When A is differential and r 
is positive and not an integer, A^ acts as a pseudodiffercntial operator (cf. [34]). 

For < T < 1, the Hilbert spaces 'D{A'^) with inner product (•, ■)r = {A'^ ■,A'^ •) 
are the interpolation spaces by quadratic interpolation between 'D{A) with inner prod- 
uct = {A-,A-) and H with inner product (•,-)o = {A^-,A° ■) (cf. [2]). This 
definition of quadratic interpolation has been used in spectral methods employing pe- 
riodic Sobolev spaces on the line in [28] . But in the finite element literature the "real 
method of interpolation" is common (cf . [12] [13] ) . Via the usual complexification of 
real Hilbert spaces both of these interpolation methods are known to be equivalent 
[IS] J [27] • Our parameter t is denoted by 1 — 6* in these references. Though inter- 
polated norms are commonly used in static {Au — /) problems and time dependent 
problems, for eigenvalue problems such methods seem only to have been used for 
rates of convergence in standard norms (cf. [7], [8]). The use of fractional powers is 
quite natural in spectral methods, and we will comment later on applicability of our 
theorems via finite element constructions. 

Now observe that by the functional calculus for selfadjoint operators, T, i.e., the 
restriction or extension of T as appropriate, is selfadjoint in 'D{A'^) for any r. We will 
use the symbol T for all of these operators, and note that in all of these spaces T has 
the same eigenvalues and eigenvectors. 

In theory the upper eigenvalues of T (lower eigenvalues of A) can be estimated 
by the Rayleigh-Ritz method in any of these spaces. The Rayleigh quotient for T in 
V{A'') is 



{A-'Tu,A''u) _ (A^ 



{A''u,A^u) (A-^u^A^u) 



for u 7^ in X'(A^). 



In practice only the Rayleigh quotients with t an integer or half integer are potentially 
useful for computation, and for r < this requires a useful formula for A~^. Our 
method is to proceed initially in the abstract, and then use duality type arguments 
to estimate convergence of standard Rayleigh-Ritz eigenvectors in other interpolation 
norms. This will resolve the convergence questions raised in the nonlinear rotating 
string example of the previous section. 
The three Rayleigh quotients 

{u,Au) ^^^^ 



{Au, Au) 



corresponding to r = 1, 



{u,u) 



a{u) {A^/^u,A^/^ 



(3.2) 



corresponding to t = 1/2, and 



^ (3.3) 

corresponding to r = 0, have been studied in the context of matrix eigenvalue prob- 
lems in 1^. Therein the quotient p.ip is termed Right Definite Lehman with zero 
shift or harmonic Ritz, l\'S.2\\ regular Rayleigh-Ritz, and (|3.3p Left Definite Lehman 
or dual harmonic Ritz. 

The Rayleigh-Ritz method in 'D{A'^), which we now denote by Ht, proceeds by 
analogy with regular Rayleigh-Ritz. One selects a family of trial vectors C Ht 

and forms the generalized matrix eigenvalue problem 



If the resulting matrix eigenvalues are denoted by 

the first monotonicity principle again guarantees that for each i — 1,2, ... ,n, Xi < 
a|"'"'. Let u^''^\ . . . , uh"'^^ be the corresponding eigenvectors orthonormalized 

in Hr-i/2- Thus Uj"''^V_i/2 = % and 

^"-)), = a("-'<5., 

where 6ij is the Kronecker symbol. 

Now let Pn,T be the orthogonal projection operator in Ht onto span{wi}f^i, the 
span of the first n trial vectors. In general, Pn,T depends on r. If {wi}°Zi complete 
in Ht, Pn,T I strongly in Ht as n ^ oo. 

Completeness of {wi}^^ in Ht is also sufficient to guarantee convergence of A^"'''^ 
to Ai as n ^ cx) for each eigenvalue of A below Aoo (cf. [37] )■ Our interest is in 
rate of convergence estimates obtainable by spectral or finite element methods from 
II (/ — P„,t) <3||r and II (J — Pn,r) Q||r-i/2i whcrc Q is a spectral projection for A onto 
the span of finitely many of the eigenvectors ui,U2, . . .. Note that while, in general, 
Pn,T depends on r, the spectral projection Q is independent of r. 

Theorem 3.1. Let Qi be the spectral projection for A onto span{ui, . . . ,Ui}. If 
{^ili^i complete in Ht, then for large n and each j = 1, . . . ,i, 

< - A, < A^^) |^|[(J-P„,00.||?-i/2 + 2A, ^E^r'j UI-Pr.,.)Qif?j 

where \\ ■ \\t and \\ ■ ||r-i/2 denote the norms of bounded operators in the corresponding 
spaces. Furthermore, there exist eigenvectors Uj^n of A in the Xj eigenspace, normal- 
ized in Ht-1/2, such that 



\\u,,n - uf-'^\l < 4A, ( 1 + ^ ) II (/ - Q|l?_i/2 + A^r^ - A„ and 



-1/2 



< 2 1 



) ||(/-P„,.)0||.- 



1/2 



where Q is the spectral projection onto the Xj eigenspace, 

S = dist{{Xj},a{A) - {Aj}), 

< K < 1, and a{A) is the spectrum of A. 

Theorem 13. II is obtained by translating the proof of the corresponding theorem in 
[33j into our context, with attention to specific constants in their estimates, and more 
detail on eigenvectors corresponding to multiple eigenvalues. The proof is included as 
the Appendix. A strength of the method is that compactness of T is not required, and 
the eigenvalue and eigenvector estimates scale correctly if A is replaced by (const.) A 
(cf. [33], [35]). But a weakness of the method is that in order to estimate Afc, 
one must employ the estimates for Ai, A2, . . . , Afc_i. Though this is consistent with 
the progression of the Rayleigh-Ritz method, in such problems as the L— shaped 
membrane considered in [1], [33| . the lowest eigenvalues have eigenvectors that are 
less smooth than those for higher eigenvalues. This gives slower rates of convergence. 
In [4], [15] localized projection estimates (though not in fractional norms), i.e., with 
Qe replaced by the spectral projection for A onto a single eigenspace of A are obtained. 
Their methods use compactness of T, and it is an open problem whether the hypothesis 
of compactness can be relaxed. 

Before proceeding to the abstract theorems for spectral implementation of Theo- 
rem l3.1l we recall a basic interpolation theorem for use in the sequel. Let Vi and Vq be 
Hilbert spaces with Vi continuously contained in Vq, written Vi Cc Vq, i.e., Vi C Vq 
with the injection from Vi into Vo continuous, and Vi dense in Vq. Then there is a 
positive definite selfadjoint operator S in Vo such that the inner product in Vi 

is given by 

{u,v)i = (Su, Sv)o, 

for all u,v E Vi — 'D{S), where (•, •)o is the inner product in Vq. For < t < 1, the 
T*'* interpolation space by quadratic interpolation between Vi and Vq is the Hilbert 
space Vr — 'D{S'^) with inner product (u, v)t = (S'^u, S'^v)o- 

Proposition 3.2. Let Wi and Wo be a second pair of Hilbert spaces with Wi (Zc 
Wo o,n-d- Wi dense in Wq . Further, let C be a continuous linear mapping of Vq into 
Wo with bound toq which is also continuous from Vi into Wi with bound mi . Then 
for each r € (0, 1), C is continuous from Vr into Wt with bound rur < m'[m]^'~^ . 

Proofs of Proposition 13.21 can be found in [5], [37]. An immediate consequence of 
this proposition is that changes between equivalent norms of each of Vi and Vo only 
affects Vr up to an equivalent norm. This follows by employing the identity map for 
both C and C'^ . 

While rate of convergence results for the linear eigenvalue problem embedded in 
the rotating string problem can be obtained by integration by parts (cf. jlOj . [llj . 
[14] ) , the fractional norm estimates needed for sharp approximation of the bifurcation 
coefficients require a more extensive theory. It is essential to be able to determine 
when a function is in the domain of a fractional power of a differential operator from 
knowledge of the smoothness and the homogeneous boundary conditions satisfied by 
the function. We now sketch results which resolve this for our purposes. 

If is a bounded domain in M" with smooth boundary results for general 
normal homogeneous boundary conditions are obtained in [2T]. Denote by H^{Vl) 
the closed subspace of H"^{Vt) consisting of the elements u S H™{rt) for which 

BjU ^0 on dn,0 < j < fc, 
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where the Bj's are differential operators of order mj < m. The boundary operators 
{Bj} are assumed to have smooth coefficients, to have distinct orders, and to be 
"normal", i.e., to have the form 

Bi = — h terms of order at most uii 

on d/dv denotes the normal derivative. Then if Vi = H^{il), Vq = L'^( fl) and 
T e (0, 1) is such that mr — 1/2 7^ rrij for any j = 1, . . . , fc 

Vr = H^^{fl) = {u G : BjU = on d^l for all j such that nij < mr — 1/2}. 

Note that for these values of r the spaces H^'^{Vl) are closed subspaces of i/™'^(ri). 
For each of the "exceptional" values r, Vr has a strictly stronger norm than that of 
H^^{Vi), the closure of H^^{n) in 

Remark. For domains with corners, let be a Lipschitzian Graph Domain in M" 
(cf. [3]). When n = 1, this means an open interval. Then for any positive integer m, if 
we let Vi be the Sobolev space H"'{VL) and Vo = = L^{n), then Vr = H"^^{n), 

< r < 1 (cf. [3] where the theory is given in the language of Bessel potentials). 
For homogeneous Dirichlet boundary conditions and Q a Lipschitzian Graph Domain, 
the following results are contained in [5D] (cf. also [H]). If, as usual, for a > 0, 
H^{n) denotes the closure of C^^n) in iJ"(f]), Vi = H^{n), and Vq = H°{n), 
Vr — H^'^{i}) whenever r S (0,1) is such that mr + 1/2 is not an integer. For 
each "exceptional" value of r the space Vr has a strictly stronger norm than that of 
H^'^{n). If, in addition, dil, the boundary of fi, is smooth enough that the selfadjoint 
elliptic operator A over defined by 

(AM,u)^2(f^) = (w,w)fl-™(n) for all v e H^{n) 

has domain _ff^'"(17) n H™{^1) (regularity at the boundary of weak solutions), then if 
Vi = H^'^in) n H^{n) and Vo = HJ^{n), Vr = H^""'{n) n H^{n) for ah T e (0, 1). 
Since 'D{A^^^) = H^{Q) (cf. [25) this resolves the interpolation problem between 
H'^"'{n)nH^{n) and L'^{fl). This statement is slightly stronger than that of Theorem 
5.3 of f2Ql, but is exactly what is proved there. As stated here, an additional possible 
application is to the homogeneous Dirichlet problem for the biharmonic operator on 
a rectangle in (cf. [22]). 

4. Spectral Implementation. Let Qq be a Hilbert space consisting of the same 
elements as Tio with an equivalent norm. We write (•, •)g and | • |o for the inner product 
and norm in Qq, respectively. Further, let M be a positive definite selfadjoint operator 
in Qq with M^^ compact. Denote by 

< /ii < < • • • < Mfe ■ • • ^ 00 

the eigenvalues of M, counted according to multiplicity, with corresponding eigen- 
vectors qi, q2, ■ ■ ■ , qk, ■ ■ ■ orthonormal in Qq. The eigenvectors {qk} will serve as trial 
vectors for the eigenpairs of A. For this purpose let i?„ be the t/o-orthogonal projection 
onto span{qi, . . . ,(/„}, and let Pn.r be, as previously, the i?7.-orthogonal projection 
onto span{qi, . . . , g„}. In the rotating string example begun in Section 2, Qq will be 
L^(a, b) with Lebesgue measure, while Hq will be the space on (a,b) with measure 
p{x)dx. 

Proposition 4.1. Assume that 

P (m1/2) Cc 2? (^1/2) (4.1) 
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and let < T < 1/2. Then if a > t and the eigenvectors Ui, M2, ■ • ■ , of A are all in 
and 

hath as n 00. Here "o" is the usual Landau symbol "little oh" . 

Proof. First note that (|4.ip implies completeness of {qk} in Then if 



\\{Pn,r-I)v\\r = \\A^ {Pn,r-I)v\\o 

= \\A^Pn,r-I) {Rn-I)v\\o 

< ||A"(i?„-/)«||o 

since Pn,T is an orthogonal projection in TCt- = 'D{A'^). Now since the identity map, 
I, is continuous from Qq into Tio and also from I?(A/^/^) into it follows from 

Proposition [3?2l that for < t < 1/2, 

||A^w|jo < c|M^t;|o, for all u £ 2?(M^). 

Herein, and in the sequel, c is a generic constant, whose value varies from place to 
place, but is always independent of the order of the approximation (n here). Thus for 



\\{Pn,r-I)v\\r<c\Kr{R^-I)v\a 

= c|(i?„-/)Ari;|o, 

since i?„ is a spectral projection for M, and hence M'^ . So since a >t and v G 'D{M'^), 

00 

\\{Pn,r-I)vrr<C ^ | (M^^^, gfe)o P 

00 

= c 5] |(M>Ar->fe)j2 

00 

00 

fe=n+l 

< c^i'^J^^" \M''v\l 0(1) as n ^ cx). 

In brief, 

II (P„,, - /)z;||2 = o (m';?") as n ^ 00. 
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Since the norm in TLr is stronger than that in Ti.r-i/2 the proposition foUows from 
Theorem O □ 
Note that we have effectively used powers of M to perform "a continuous abstract 
integration by parts". Proposition 14. II provides potentially useful information for r = 
or r = 1/2. But for < r < 1/2, the Rayleigh quotient is usually not computable 
since typically only complicated expressions for equivalent norms arc known, and that 
at best. 

We now give the analogue of Proposition 14. II for the harmonic Ritz method, i.e., 
the Rayleigh quotient with t — 1. Qo and M are as previously defined. 
Proposition 4.2. Assume that 



V{M)Cc'D{A) (4.2) 
and let < T < 1. Then if a > t and the eigenvectors Ui,U2, ■ ■ ■ ,U£ of A are all in 

A;"-) = A,+o(Mr+?^),l<j<£, 

and 

\\Uj,n~U^;^^^\\r^o{fi:-l), 1<J<£, 

both as n —^ oo. 

The proof directly mimics that of Proposition 14.11 As there, potentially useful 
information has been obtained for r = 0, 1/2, and 1, and for t — 1, this is applicable to 
the nonlinear string problem. But in that problem, the critical value of r to treat the 
expansion of the nonlinearity is r = 3/4. The following theorem enables us to apply 
regular Rayleigh-Ritz, i.e., r — 1/2, with t large enough to treat the nonlinearity. 
This requires stronger hypotheses that are only satisfied when A has compact inverse. 
An example in another context [5] applies hypotheses like those of the preceding two 
propositions to A without compact inverse. 

Theorem 4.3. Let t > 1/2, assume that 'D{An) = 'D{A'^), with equivalent 
norms, and further assume that 'D{M^'^~^/^) Cc 'D{A^'^^^/^). Then if for each j = 
!,...,£, both 

Eoo 4r — 1 I / \ 12 

Mfe \{Uj,qk)o\ ^/ 2r-n ,.o\ 

and Uj e T>(M'^) with a >2t — 1/2, one has 

11 (n,l/2)i| / T — a\ ■ 1 n 

\\Uj-u) ||r=o(/^„_^J as n ^ oo,j ^ 1, ... ,t. 



Proof. Let uj""* = m^"'^^^^ and define by 

/ll-2r (") 

i.e., 

u;,,„ = A^^-i (u, - uf'^) . 
11 



Our hypotheses guarantee that wj^n is well defined, and that the following manipula- 
tions are valid. Then 

i (") \ 

The Rayleigh-Ritz method is such that if = for some n, it vanishes for all 
greater n, and there is nothing to prove. Otherwise, by the definition of Wj^m 

Noting that 

\\w3,n\\l/2 ^ \\A^~^''Wj.n\\2T-l/2 

it remains to prove that the latter ratio is 0(/i^_|_i^^), as n —> oo, since by Proposition 
14.11 \\uj — Mj"^||i/2 = as n — * CX3. For this purpose observe that 

since V{A^) = X'(M^) and ^(A^^-i/^) Cc V{AP^~'^/^). Thus it suffices to show that 



|M2— 1/2„|2_ \M^-^/^x\l 



\Mr 12 ^ ™P I 12 

= 0{fl^'^^^) as 71 OO, 

where we have used the fact that the q^s are eigenvectors of AI, and hence of M"^. 
Then if R„ is the spectral projection for M onto span{qi, . . . , g„}, 

sup --2 

xGspan{AI'^ Uj ,qi ,....qn} I'^lo 

sup - 

^ |M--V2j;„^|2 \Ar-y^i-R^)x\l 

~ 1 I 12 ' I 12 

xespan{(I-R„)M^Uj,qi,...,q„} \ F|o F|o 

by the Pythagorean theorem. The latter is dominated by 

|M-i/2i?„x|2 |M--i/2(/^i?„)a;|2 

sup j— 12 + sup j— 12 

} Flo 

which is equal to 

|M-i/2j?„^|2 ^ |M-^V2 (/_ j^JM-^.,■|2 

xespan{qi,...,qn} li^ ^ Rn)M'^Uj\Q 

12 



since i?„ and / — i?„ are complementary spectral projections for M , and therefore 
for the powers of M. Since r — 1/2 > the remaining sup is, by the variational 
characterization of eigenvalues of M'^^^^'^, exactly 

|M--V2 5„|2 



2 — H-n ■ 



The latter ratio is 



Eoo 4r — 1 I / \ \7 

fc=n+l Mfc \{uj,qk)o\ 

which by hypothesis is 0{fi'^^^^). Since fin+i > A^n and 2t — 1 > 0, this concludes the 
proof. □ 

Note that use of regular Ray leigh- Ritz, rather than "r-Ritz" with r > 1/2, does 
not change the "optimal" rate of eigenvector convergence in the r-norm as given in 
Proposition 14.21 Since, when applicable, dual harmonic Ritz gives faster eigenvalue 
rates of convergence for a given set of trial vectors we now present the corresponding 
theorem for this case. 

Theorem 4.4. Let r > 0, assume that 'D{M'^) = V^A^), with equivalent norms, 
and further assume that T>{M'^'^) Cc 'D(A^^). Then if for each j — !,...,£, both 



and Uj £ 'D{M'^) with <t > t, one has 



0(m„+i) asn^co, (4.4) 



(n,0)|| f T-(j\ 

u) '||^ = o(^„^J as CO, 3 



The proof of Theorem 14 . 41 mimics that of Theorem l4.3l starting by defining ujj -yi as 

A~'^'^ Wj^n = The corresponding theorem using a harmonic Ritz calculation 

with r > 1 would entail defining Wj^n by A^^'^'^ Wj^n — Uj — We will not pursue 

this. 

We now examine convergence of Rayleigh-Ritz eigenvectors in the Hq norm, 
weaker than the energy norm. The following is an abstraction of a duality argument 
from finite element theory as presented in [13j . 

Theorem 4.5. Assume that 'D(A) = 'D{M) with equivalent norms |Mu|o and 
\\Av\\o. Then if a > 1/2 and v G V{Ar), 

\\v - Pn,i/2v\\o = oiUnli) as n-> CO. 
Moreover, the conclusion of Proposition \4-l\ can be supplemented by 
\\uj - uj"'^/^^||o = o(/i^^i), l<j<e, as CO, 

provided that the first £ eigenvectors, Ui, . . . ,ui are in T>[M'^). 

Proof. For simplicity of notation denote Pn,i/2 by Pm and for v G T>{M'^) let w„ 
be the unique solution of 

AWn =V - PnV. 
13 



If Wn = there is nothing to prove. Otherwise, 

\\v - PnV\\l = (w - PnV,AWn)o 
= a{v - PnV,Wn) 
= a{v - PnV, Wn ~ PnW„). 

So by the definition of Wn, 

\\v - PnV\\l/2 \\Wn ~ P„W„||i/2 



\v - PnV\\o < 



\AwJ\o 



\\V - PnV\\i/2\\Wn - PnWn\\l/2 

\U \ • 

\MWn\o 

Since Wn G T^iA) = 'D{M)^ it follows as in the proof of Proposition 14. II that 

||W„ - PnWn\\\/2 < C | | q , 

and since v £ 'D{M'^) with a > 1/2, it follows similarly that 

11^^ - Pnv\\l/2 < cfil+l'' \M-v\l 0(1) as n ^ OO. 

So (1131) yields 

\\v - PnvWo < cfi-l, fi^-l^ \M''v\l 0(1) = oifi-l'[) as n ^ oo 

and the theorem follows from Theorem 13. II □ 
It is possible to use Proposition 13.21 to obtain \\v — Pn.i/2^\\T = o(^^j^i) with 
0<r<l/2asn^oo under the hypotheses of Theorem 14.51 But this does not 
immediately extend to a Rayleigh-Ritz eigenvector estimate in the r— norm. That 
remains an open problem. The theorem for harmonic Ritz eigenvector estimation in 
the energy norm is as follows. 

Theorem 4.6. Assume that'D{A^/'^) — P(M^/^) with equivalent norms \M^/'^v\q 

ana 

11^3/2 

v\\q. Then if a > 1 and v £ 'D{M'^), 

\\v - Pnsv\\l/2 = oifil-ll") asn-^oo. 
Moreover, the conclusion of Proposition \4-^ can be supplemented by 
\\uj - Uj"'^^||?/2 = oi^il-f^"), l<j<£, as n oo, 

provided that the first £ eigenvectors, ui, . . . ,ue are in 'D{M''). 

The proof of Theorem 14.61 mimics that of Theorem 14.5) provided one starts by 
defining w„ by A^^'^ Wn = v — Pn,iv, and so is omitted. It is then possible to use 
Proposition 13.21 to obtain — P„.if ||t = o(/iJ^^'[) with 1/2 < r < 1 as n — s- oo, but 
again this docs not lead directly to a corresponding eigenvector estimate. 

5. The Model Nonlinear Problem, Part II. There are at least two ways to 
formulate the spaces Hq and Qq, and the operators A and M to implement the theory 
of Section HI for this problem. It is advantageous for our purposes to write (|2.5p as 

Avo = -^Vq ^ i^ovo, vo{a) ^ vo{b) = (5.1) 
14 



with V{A) = H'^{a,b) n H^{a,b). Then A is selfadjoint in L^{a,b) with measure 
p{x)dx, which we take as Tig- The quadratic forms which are needed for harmonic 
Ritz and regular Rayleigh-Ritz are then found to be 

{Av,Av) = / ^^dx on V{A), 

Ja P 



{Av,v) ^ / {v'f dx on V{A^^^) = H^{a,b), 



and 



As AI we take 



b 

2 



v,v) = V pdx on TYo- 



Mv^-v", v{a)^v{b)^0, (5.2) 

selfadjoint in Qq = L^{a,b) with Lebesgue measure, dx. Then = T>{M) with 

equivalent norms, and the eigenvectors of M generate the usual Fourier sine series. 
Note that Av = at a and b if and only if v"{a) = v"{b) — 0, and since the same is 
true for M, we have 'D{A^) = V{M'^) (with equivalent norms.) But if w e 'D{A^), 

A^v = -- (--v' 
P \ P 

=K^"""'G)'°"'"G)""") 

must vanish at a and 6, while if u e I?(Af'^) 

must vanish at a and b. Independence of the boundary functionals shows that 'D{A'^) 
and 'D{M^) are different closed spaces of H^{a,b), unless p'{a) = p'{b) = 0. These 
boundary functionals are unstable in the topology of iJ^/^ (a, 6) (cf. [H])- So quadratic 
interpolation between ViA^) and 'D{A^), and between I?(Af3) and 'D{AP) (as sketched 
in Section [3]) gives 

ViA") ^V{A'r) for aU < 2 + i ^, 

with equivalent norms. So one can apply Theorem 14.31 with 2r — 1/2 < 9/4, i.e., 
r < 11/8 to our problem. 

The other standard way to treat (|2.5p is to let zq — ^ to get 

Azq^ — ^ ( ^ ) ^vqZq, zo(a) = zo(6) = 0. 

The so defined operator A is selfadjoint in L^{a^ b) with the usual Lebesgue measure. 
So taking this space as both TYq and Qq, and M as above yields the framework of 
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Section 21 But an analysis like that in the previous paragraph yields only 'D{A'^) = 
'D{M'^) for all < 1 + | = |- This yields slower rates of convergence, though Theorem 
Hillcan be applied with 2r - 1/2 < 5/4, i.e., r < 7/8. 

We now proceed to apply the convergence theory of sections [3] and [4] to solution 
of (|2.5p and (|2.6p with the operators A and M chosen according to (|5.ip and (|5.2p 
respectively. As is well known, the eigenvalues (all simple) of M are given by 

^J-k — /, NO , K — 1, 2, ... , 

(o — ay 

with corresponding eigenfunctions (unnormalized) 

kT:{x — a) 
qk = sm — ; , a < X < b. 

b ~ a 

Since fik is proportional to fc^ and 'D{A^) = V^M^) for all < 9/4, Proposition 
14.11 shows that the Regular Rayleigh-Ritz convergence rates for any index j are 

^(n,i/2) ^ _^ o(n"") for aU a < 7, 

and 



(n,l/2) 



II 1/2 = o(77r") for all a < 7/2. 



But the result we need for calculation of the bifurcation coefficients is that of Theorem 



In order to accomplish this we will need to verify that (|4.3p holds. This will be 
carried out via integration by parts, which will also slightly improve and simplify 
the statements of the estimates of the preceding paragraph. Our abstract "small 
oh" estimates are analogues of a classical theorem on Fourier series coefficients of 
smoothly periodic functions, while the improvement comes from the corresponding 
"big oh" estimates for functions satisfying Dirichlet conditions (cf. [32], p. 130). This 
use of integration by parts is also employed in TO], [11], [T4|. While the abstract 
theory of Section [5| conveniently identifies the possible "big oh" power, it does not 
directly imply such improvement, as follows from examples in a different context ([20], 
pp. 155-156). 

So let (A, u) be an eigenpair for the selfadjoint operator A, and observe that from 
the proof of Proposition |42| that for < t < 1, 



||(p„,,-/)u||,<c|(i?„-/)Aru|o. 

{Rn — I) M'^u is just the tail of the Fourier sine series of M'^u, 



oo 

\ ^ 2t 



|(i?„-/)Aru|o= 2^ \{M-u,qk)oV^ 2^ ^^k\iu,qk)or. 

k—n+l k—n+1 

This is also the quantity of concern for larger values of r for application of Theorem 

Since normalization of qk in Qq = L^{a, b) is independent of k, we can bury it in 
the constant c, and integrate by parts six times to obtain 

kii^x — a) 



u sm ■ 



■ dx 



- f (.»(6)cosfc7r-.™(a)- sin " dx] . (5.4) 
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The boundary terms in the other five integrations by parts all vanish. To see if this 
predicts the rate of convergence, first note that since the eigenvector u is in the domain 
of any power of A, (|5.3p implies that 

= -2p (-^ u'" at a and h. 



Then since [ -u"\ — —\u' '\t follows that 



( — ) u' at a and h. 

P 



Now, u' {a) 7^ 7^ ^'(^) since u is an eigenfunction and u{a) = u{b) = 0, by uniqueness 
for the initial value problem. Thus w™(a) 7^ 7^ unless p'{a) = p'{h) = in 

which case one obtains a more rapid rate of convergence. Excluding this special case, 
we have for < t < 1, 

/>OC 

k=n+l 

since pk — (^f^Ja)^ ■ ^'^ follows as in the proof of Proposition 14.11 that 

xi^^m = x + o{n-'), 

and 

and from that of Theorem 14.51 that 

all as n — > cx). In addition, (j5.4p implies that there exist constants ci and C2, inde- 
pendent of k, such that 

< \{u,qk)a\ < C2k-^ 
for large k. So the numerator in (|4.3p is dominated by a constant times 

00 ^8t-11 



fc=r! + l n+1 



n 



11 - 8t 



provided that 8t— 12 < —1, i.e., that r < 11/8. Recalling that we must have a < 9/4, 
this is consistent with the restriction r < § + -j in the statement of Theorem 



n+l 9 - 4t 



Similarly, the denominator in (j4.3p is minorized by 

00 00 „oo -I- n4r-9 

^ A;4r^-10^ ^ / ^4r-10^^ ^ + 

and so the estimate (|4.3p is 



sr-11 \ 
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as n — > oo, for 1/2 < T < 11/8. 

Since we have already verified that \\u — w'-"'"'^^^-' |1 1/2 — 0{n~'^/'^), it now follows 
from the proof of Theorem 14.31 that 



^(n,i/2) = A + 0(n-^), and 

(5.5) 

for 1/2 < T < 11/8 as n ^ 00. 

We now examine use of Rayleigh-Ritz approximate eigenvalues and eigenvectors 
to calculate vi as given in (|2.8p . Therein, [vq, vq) is an eigenpair of Au = \u. We will 
write (A'"',it("^) for the corresponding Rayleigh-Ritz eigenpair (A'"'^/^^ u*^"^^/^'). 
Thus we estimate the rate of convergence of 

A(") (^p{x) , to A {p{x) u {u'f\u) . (5.6) 

Note that since we already have convergence rate estimates for A'") to A, it suffices 
to estimate the difference between 

/ p(a;)w(")Vu(")') dxand / p{x) {x) {u' {x)f' dx . 

Write the difference between these integrals as 

j\{x)u^''''>^ (^(u^-^^y' ~{u\x)f'^dx + p{x) (u'(x))'* u^) ds. (5.7) 

Since u' is smooth, the second integral in (|5.7p is dominated by 

sup \p{x) {u'{x)f\ / -u'\dx, 

x^{a,b) Ja 



and by the Cauchy-Schwarz inequality the latter integral is dominated by 



So by the preceding eigenvector estimates the second integral in (|5.7p is 0{n ^/^). 
Since u'"^ is uniformly bounded, the first integral in (|5.7p is dominated by 



sup 



|p(:r)u(")'| / {{u^^^yy ~{u'{x)f\dx, 



and by the Cauchy-Schwarz inequality the latter integral is dominated by 

/ b \ I I) 1t — l 

IJ^ |«(")' - ufdx\ /j E ("'"^T (^')'*"-'"' l''^^ 
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Fig. 5.1. Rate of convergence of \{n) to the first eigenvalue X = tt^ + 1/4, The equation of the 
linear fit is given by y = 4.4156 - 7.22504x. Note that |A - A(20)| = 2.795 X 10"*. 
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Fig. 5.2. Rate of convergence of the bifurcation coefficient to ui for t = 1. The equation 
of the linear fit is given by y = 8.14125 - 7.29132x. Note that \ui - | = 1.86 X lO-**. 



19 



t 




(20) 


1 




_4(e-l)7r''(7r2 + i)(7+47r2) 


-18.008997020330582 




e(l+20Tr^+647r*) 


2 




-l)7r''(7r^ + i)(437+8247r^ + 1447r^) 


-75.15014087855786 




64e^(l+147r^+497r*+367rB) 


3 


-16(e^-l)7r**(7r^ + i)(1709+55407r^+3856ir*+320ir'') 


-571.7347727528597 


96^^(7294 


97207r^ +393127r«'+524807re +163847rS ) 



Table 5.1 

Exact and calculated values of vi for i = 1, 2, and 3. 
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Fig. 5.3. Rate of convergence of the bifurcation coefficient u\ 
of the linear fit is given by y = 10.8453 — 6.97158a;. Note that \v\ 



(n) 



to u\ for 4 = 2. The equation 
4.3 X 10-^5. 



(20), 



By the Sobolev inequality |22j . (15. 5|) with any r > 3/4 imphes uniform conver- 
ffence ofw^") to u' , so the latter factor is bounded. Then the eigenvector estimate 
(j5.5p with r = 1/2 gives a rigorous rate of convergence estimate for the first factor of 
0(n^^/^). In other words, the rate of convergence estimate for the Rayleigh-Ritz ap- 
proximation to the first bifurcation coefficient is the same as that for the eigenvector in 
the energy norm, for any t — 1,2,.... The following calculations seem to indicate that 
our convergence rate estimates for the first bifurcation coefficient are quite conserva- 
tive. This may be due to the need to employ the Cauchy-Schwarz inequality in the 
absence of estimates (see Fig. 5.1-5.4). For the problem at hand, with p{x) = l/x^, 
a = 1, and b — e, the exact and calculated values of i^i for t = 1,2, and 3 are shown 
in Table 5.1. The eigenvalue and normalized eigenfunction are Ai = tt'^ -|- 1/4 and 
Vl (x) = sin (In x) . Calculations were performed using Wolfram Mathematica 6. 0. 
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Fig. 5.4. Rate of convergence of the bifurcation coefficient i/^"' to ui for t = 3. The equation 
of the linear fit is given by y = 12.4435 - 6.26615x. Note that \ui - u[^°^ | = 1.1 X IQ-^. 



We now investigate the Rayleigh-Ritz approximation for vi. Convergence of our 
procedure for i/2,V2,--- will then be clear. While this will verify convergence, the 
method about to be presented appears to be impractical as a numerical algorithm for 
these higher order terms. 

Recall that by (12. 6p vi is the solution of 



v'l - uq p{x)vi = vi p{x) «o + ^ p{x) vq (wq)^* , vi (a) = vi [h) = 0, 



where by ([2 



with (wq, vq) — 1. With this value of vi and normalization of vq, this is equivalent to 



or 



where A — vq has been inverted on the Tip orthogonal complement of vq, and the 
component of vi parallel to vq is taken to be zero. We now fix attention on the least 
eigenvalue of A, i.e., vq = Ai and vq = wi, the procedure for higher eigenvalues being 
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entirely analogous. The spectral expansion of (|5.8p is then 

i=2 

N2t 



V"';" -"''-u,. ,5.9) 

Recalling the Rayleigh-Ritz normalization 
where we again write w^"^ for etc., gives 

, ^ / iri) ( (n)\ 

(„) _ Aj_ ^ \ ^ V ^ J 'Wo „(n) 

3=2 

As noted previously, by Theorem 14.31 Uj in Co[a,6] for each j. Thus we 

have term by term convergence in Ti.i/2, and convergence of wj"^ to vi in Ti.i/2 '^iH 



follow from the dominated convergence theorem for vector valued functions (cf. [17] ) 
provided that ||wJ"''||i/2 is uniformly 
Rayleigh-Ritz normalization we have 



provided that ||wJ"''||i/2 is uniformly bounded with respect to n. Again using the 



This implies that 




'"•(»s"'T'."r>o 



A<"> 



2 / \ 

A") 



^ II / /\ 2t 1 1 2 

< c||ui [ui ) llo, 

and convergence in 7ii/2 is proved. Rates of convergence in various interpolation 
norms for a fixed finite number of terms in uJ"-' to those in vi follow from the preceding 

(n) 

theory. But this is not sufficient for estimates of the rates of convergence of Vi to 
vi. Theoretical convergence of 1^2"^ to 1^2, ^^2"^ ^'^^ ^^'^^ follow similarly. 

6. Examples and Remarks. Let V denote the unit ball on M'' and consider 
the branching problem 

— Ay + yP — X p[x) y in P, y = on 91?, 

where A is the Laplacian, p is an integer greater than or equal to 2, and the density 
function p{x) is smooth and strictly positive on T>. As usual dV denotes the boundary 
of V, \x\ = 1, and V is the closure of V, \x\ < 1. By setting y — u, with e small 
and positive, we seek normalized solutions to 

- Au + euP = Xp{x)u inV, u = OoiidV. (6.1) 
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A variant of (|6.ip for a bounded simply connected domain fl in M'^ is obtained for 



A u 



fl, u on dfl, 



by use of a conformal map w ~ f{z), U{z) = u{w), to the imit disk V obtaining 



-AU + e 



dw 


2 


dw 


2 


iJP = A 


U inV, 


dz 




dz 



= on dV 



(6.2) 



(cf. [25]). Provided that is smooth so that / is conformal on dfl, p{x) ~ \^\^ 
has a positive lower bound and our theory applies to ()6.2|) as well as to ()6.1|) (with 
the exception of strictly radial solutions of (16. ip ). 

Existence of solutions of (16. ip or (|6.2p bifurcating from a simple eigenvalue of the 
linearized (e = 0) problem follows as in [31 by use of the Banach space C^'"(2?) of 

functions having Holder continuous derivatives with exponent a, < a < 1. The 
formalism sketched in the rotating string problem leads to the bifurcation coefficient 
for (O) 



iP+i dx 



(6.3) 



where u is the eigenvector corresponding to the simple eigenvalue A of (j6.ip with e set 
equal to zero and normalized by J^, p(x) dx — 1. An analogous formula holds for 
(|6.2p . and we assume that p{x) is such that u is unknown so that we must approximate 
u appropriately to approximate i^i. 

The interpolation theory is directly analogous to that previously detailed in the 
rotating string problem. The eigenvalue problem for the operator A is given by 



Avq = Avq in V, vq = on dV, 

P 

with P(A) = H'^{'D) n HliV). The operator A is selfadjoint in L'^iV) with measure 
p{x)dx. As AI we take 



Mv 



-Av in 



on dV, 



which is selfadjoint in Qo = L?'{T>) with Lebesgue measure, dx. As previously, 'D{A) = 
I?(M) with equivalent norms, and since Av = on dV if and only if —Aw — Mv — 
on &D, we have 'D{A^) = Vi^A-P) with equivalent norms (cf. [571 for the regularity 
theory). But if u e V{A^) 



-'-A 
P 

1 /I 



1a 

p 



■ A^v + 2V 



V{Av) 



A ( i I Aw 



P \P \P 
where V = gradient, must vanish on dV. On dT), the latter is 



and if t) e V{A'P), 



A^v = - - A% + 2V - • V (Av 
P \P \P. 



M^v = A^v 
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(6.4) 



(6.5) 



must vanish on dV. Since v E T^iA) = 'D{M), Av = on dV, and so V (Au) is 



orthogonal to dV, and (|6.4p and (|6.5p are the same boundary conditions if V yj^j is 

tangential to dT> (or vanishes there). We now show that when this is not the case 
the boundary conditions (|6.4[) and (|6.5p are not equivalent. This leads to the same 
critical value r = 9/4 as in the rotating string example. 

So suppose there is a point P on dD where ^ (^p) ' ^ 7^ 0, where it is the unit 
outward normal to dT>. By continuity there is a relative neighborhood M oi P in dV 
on which V • 7? is non-zero of a fixed sign. Now note that if u is an eigenfunction 
of A with eigenvalue A, the boundary condition A^u = becomes 

A^u = 2\pV (-] ■\/u = V (-] -n^ ondV, 
\Pj \Pj on 

since Vu is orthogonal to dV. If there is a point of M at which 7^ there is a 
relative neighborhood Af' C A/" of P on which ^ is non-zero of fixed sign. Then the 
boundary conditions (|6.4[) and (|6.5p define different closed subspaces of H^iV), i.e., 
'D{A^) ^ V{M'^). The other possibihty is that |ii = on AA. But then u has zero 
Cauchy data on M . Since T) is strictly convex at P, by [3^ there is a disk 5 centered 
at P with u = on P n 5. But then by a theorem of [55] (cf. also [30]), u = on I?, 
contrary to u being an eigenfunction of A. 

So if V is not tangential (or zero) at any point of dV the results of [3T] give, 

by interpolation between V{A?) and V{£-), and between Vi^W^) and ^(M^), that 
X'(yl'') = X'(A/'') for aU < 9/4. By details in [H] not elaborated on in Section H 
V cannot be increased to 9/4. Direct application of Proposition 4.1, Proposition 4.2, 
Theorem 4.5, and Theorem 4.6 now proceed exactly as in the rotating string example. 
For example. Proposition 4.1 gives the Rayleigh-Ritz eigenvalue approximation 



A^"''/') = A, + o(M^;f ) for all u < 9/4 
as n — !■ 00. However, in R'', /i„+i ~ cr?!'^ (cf. [38J) so that 

aJ"'^/') - Aj + 0(^2(1-2-)/'') for aU a < 9/4 

as ri — > 00. 

Before looking at rates of convergence for Rayleigh-Ritz or harmonic Ritz approx- 
imations to the bifurcation coefficient (|6.3p we look, as in the rotating string example, 
to improve these estimates to "big oh at 9/4" by multi-dimensional integration by 
parts. With u an eigenfunction of A, i.e., 

— Au — A p(x) u in 2?, u — on dV, 
and qk an eigenfunction of M, i.e., 

- Aqk = ^^k Qk in V, qk = on dV, 
we recall the Lagrange-Green identity 

J (uAqk - qkAu) ^ J ~ ^''^) ^'^'^'^ 
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Herein, as previously, ^ is differentiation in the outward normal direction to dT). By 
first substituting A^u for u, and then Au for u in (j6.6p . we obtain 

uqkdx — — [ Ik A^udx H — g / A^udS. (6-7) 

To obtain (|6.7p we have also used u = Am = (jf^ = on dV. (|6.7p is the multi- 
dimensional analogue of (5.4), and since = Vg^: • ~n and div(gradqfc) = Ag^., we 
expect that 

uqkdx ~ 0{fj,f.^^^) as k ^ oo, 

as in the one dimensional case. But verification of (4.3) is problematic, unless we 
restrict to the radially symmetric case, i.e., that p depends only on the radial variable. 
The problem arises with the denominator of (4.3). 

Now we treat the two dimensional case in some detail, with later comments on 
higher dimensions. In the usual polar coordinates, the eigenvalue problem for qk is 
now 

d^qk . I dqk , 1 d'^qk\ r ^ . r, r , 

fik qk for r < 1, = tor r = 1. 



g^2 J, Qj, J.2 QQ 

Separation of variables yields the eigenvalues 

P'm,s — jm,sJ 

where jm,s is the s*'* positive zero of the Bessel function J^. These must be ordered 
in an increasing sequence to obtain the /ifc's. For m — 0, the eigenvalues Jq ^ of M 
are simple with normalized radial eigenfunctions 

qo.s{r,0) ^ ^ -rj-. — ^, s = l,2,..., 

while for m — 1,2,... the eigenvalues j, of M are double with normalized eigen- 
functions 

'LU'T. &)^\ - -J r — V ^ 

V TT Jm+l(jm.,s) 



and 



2 Jjnijm.s'^^ . , 

— sm m f 

Jm-\-iyjm.s) 



for s = 1,2,..., (cf. [T], |23|). Since the outward normal derivative is given by 
differentiation with respect to r at r = 1, we get 

dqo,s I (-1)-+^ . 19 



dr ^ 
and for m = 1,2,..., 



dr ''■-1 V TT 
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jm,s (-1)'+^ COS m 61 



and 



for s — 1,2, . . .. So the expectation that (|6.7p gives uq^dx — 0{pL^'^^'^) is borne 
out. 

In the radial case we have = jo,s = {s — \)t^ + 0{^) as s — ^ oo, fT] and an 
analysis like that in the rotating string example shows that (4.3) is satisfied with 
i < T < So in the radial case, using the fact that the Sobolev inequalities (cf. 
[22]) imply that convergence in Tir for any r > 1/2 implies uniform convergence. 
Theorem 4.3 and Theorem 4.5 give convergence of the Rayleigh-Ritz approximation 
to vi in (6.3) at the rate of 0(n^^/^). This is the same rate of convergence as that 
for the Rayleigh-Ritz approximation to the eigenvector u in the norm. 

In the non-radial case one can invoke Proposition 4.2 and Theorem 4.6 to obtain 
the convergence rate 0(/i„^(*) = 0(n^^/*) for the harmonic Ritz approximation to 
vi, since the eigenvector approximation is uniform. The Rayleigh-Ritz method can 
also be employed via the Sobolev inequality 

H\V) Cc W^iV) Ce W°{V) = L^iV), (6.8) 

where q > 2 and t = 2/q (cf. [22]). So convergence of u*^"^ to u in H^(T>) implies that 
is bounded in L'^iV) for any q > 2. Thus Holder's inequality yields convergence 
of the Rayleigh-Ritz approximation to vi in i«'(X»), i ^ = 1. This is majorized by 

convergence in L^iV) where the rate is at least 0{fi~^^^) = 0{n'~^^'^) as n — > oo, by 
Theorem 4.5. 

For dimensions higher than two one may employ spherical harmonics and Rayleigh- 
Ritz as above for p — 2 and p = 3 in dimension three, and for p = 2 in dimension 
4. The limitations on p stem from the need for t > in the analogues of (|6.8[) [22] . 
Harmonic Ritz is directly applicable for any power p in dimension three via uniform 
convergence. Further conclusions can be obtained from Theorem 4.3 and Theorem 
4.5 in the radial case. One can also consider "derivative nonlinearities" such as |Vup 
in place of in (|6.ip for further examples, and the addition of linear potential terms 
q{x)u. The dimensional dependence of the various conclusions comes from the Sobolev 
inequalities. 

We now close with some remarks. For domains with shapes other than balls the 
main constraint on the preceding methods is the availability of computable approx- 
imate eigenvectors, the g^'s. Such are obviously available for rectangles, or for two 
dimensional domains mapped conformally onto rectangles. But the interpolation re- 
sults of [21] are not application. Theorems of [20] do apply in the case of homogeneous 
Dirichlet boundary conditions for A, but only up to t = 1. So the theory of quadratic 
interpolation between closed subspaces of Sobolev spaces determined by homogeneous 
boundary conditions is currently inadequate to obtain results corresponding to the 
previous over domains with corners. 

The finite element method is of course well adapted to the construction of approx- 
imate eigenvectors over domains with irregular shapes. One estimates the projections 
in Theorem 3.1 by use of the "approximation theorem", as detailed for r = 1/2 in 
|33j . Finite element versions of Theorem 3.1, Proposition 4.1, Proposition 4.2, The- 
orem 4.5, and Theorem 4.6 are attainable. But Theorem 4.3 and Theorem 4.4 each 
depend on analogue of the "inverse inequality" of the spectral Galerkin literature (cf. 
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|28|). It is an open problem as to whether the inverse estimates of finite element the- 
ory [13] can be similarly applied. Quadratic interpolation between closed subspaces 
of Sobolev spaces determined by homogeneous boundary conditions in the context of 
the finite element method is considered in 12J, whose results supplement those of [20] . 
and in 43j. In the latter reference the effect of the boundary conditions in quadratic 
interpolation is not presented, but it does not affect the rates of convergence due 
to the construction of the finite elements. It does, however, affect the norms of the 
interpolation spaces in the "exceptional" cases. 

We have concentrated on the problem of estimation of bifurcation coefficients. But 
the theorems of Section 4 are also germane to estimating eigenvalue bounds by the 
so-called "eigenvector free method" , EVF for short. These bounds are complementary 
to the Rayleigh-Ritz bounds, and hence provide explicit error estimates. The EVF 
technique was originated in [18j , and our results herein relate to the convergence 
theory of [5]. This problem remains to be explored. 

Appendix. We give the proof of Theorem 3.1 in the case t = 5, i.e., regular 
Rayleigh-Ritz. This enables us to simpHfy notation by writing || ||a for || \\i/2, 

Pn for P„_i/2, and A^"'* for A^"'^''^\ The theorem as stated is obtained simply by 
substituting r for 1/2 in the appropriate places. Note that for r < 1/2, this replaces 
Ti = Ho by 7i = T-It-i/2, a "negative norm" space. 

Proof. Let £1 = span{ui, . . . ,ui} and let ee be the set of vectors in £{; which are 
normalized in Ti. The first step is to prove that if 



.(") 



then 



— max \ 2Re{u, u — Pnu) — ||u — P„w||^ I < 1, 



(jA.ip will follow easily from the minimax principle provided that Pn£i is ^-dimensional. 
To see this, observe that if Pn£i has dimension less than there exists w ^ et with 
PnW = 0. But then 



4"^ > \2Re{w,w-Pnw) - \\w~Pawf\ = \\wf = 1. 

n large 



(n) 

Thus, if < 1, which by completeness of the Rayleigh-Ritz trial vectors is true for 



.(n), a{v) ajPnu) 

A) < max -j— j-r = max—- — —. 

Since P„ is an orthogonal projection in I?(a), a{P„u) < a(u), and also for u in 

||P„U||2 \\U-{U-P„U)W' 

= \\u\\'^ — 2Re{u, u — Pfiu) + ||u — P„u|p 



Thus if p1"^ < 1, 



A(") < max = ^ 
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or equivalently, 



the corresponding estimate for eigenvalues of T. 

The second term in the definition of p^"-* is dominated by || (/ — P„) Q^jp, where 
is the orthogonal projection onto S^. To obtain an analogous estimate for {u,u — 
Pnu), first write u £ eg as u = Y^l^i CiUi. Then 

{Ui, U - PnU) = \Y^a(ui,u - PnU) 

= XY^a{ui - PnUi, u - PnU), 

since P„ is an orthogonal projection in X'(a). Thus if w G 

e 

{U, U - PnU) = ^ CiX~'^a{Ui - PnUi, U - PnU), 
i=l 

and so, 

e 

2\Re{u,u- Pnu) \ = 2 1 ^ ReciX^^a{ui - PnUi,u- P„u)| 

i=l 

t 

< 2\\U - PnUWa WCiXr^ {I - Pn) «i||a 
i=l 

Now, since u G eg, 

\\U-Pnu\\a = II (/-P„)Q^M||a 

< ||(/-Pn)Q^||a|h||a, 

e e 



||l^-PnW||a< A^'^'||(/-P„)Q^|U. 



^ ||C,A-1 (/ - Pn) U,\\a < J2 l^^l V'll U - ^n) g^tiilU 

<||(/-P„)g,|U^|c,|Ari||«,||, 

i=l 

= 11 iI-Pn)Qe\\aJ2\^^\X;'^' 
i=l 
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and so. 



Finally, 



since — A^^^^. The Cauchy-Schwarz inequality now yields 

max2|i?e(u,ii-P„u)| < 2Ay' V A^^/^ || (j _ p^^) g^||2 

1=1 

since l^jP = 1- Thus we have the eigenvalue estimates 

< - A, < {^{1 - PMf -r1\,(^\T^ W-PM\t\^ 
or equivalently 

0<^-^<^ |^ll(/-Pn)g.|P + 2A,(^^Ari^ ll(/-i'n)0.||^j . 

To derive the eigenvector estimates, first note that — ||'«^"''|| = 1, so 

aiui — w^"'') — aiui) — 2Re a{ui, u^"^) + a(M^"'') 
= Xi-2XiRe (M,,4"^)+Af' 
= A^ (2 - 2Re {u,, 4"^)) + ^ - A^ 

= A, - 2i?e (zi,, 4")) + + A^") - A, 

= A,h,-4")|P + A("'-A,. 

So to estimate the error in energy it remains to obtain an estimate for 
\\ui — To do this we will need the simple identity 

(a^") - A,) = a{PnW,uf) - Xe{BaW,uf) 

= A^(ii^-P„ii^4"^), (A.2) 

(n) 

valid for any normalized eigenvectors ug corresponding to Xg, and uj, corresponding 
to A^"^ For use in the sequel, please observe that in (|A.2p the eigenvector ue can be 
selected in an n-dependent fashion in the Xi eigenspace. 

So let A^ — A^+i = . . . = Xi+K have multiplicity K + 1, and define (5 > by 

S^dist{{Xi},a{A)-{Xi}). 

Choose < K < 1. Then by the preceding there exists N, depending on k, such that 
for all n > iV and j such that Xj 7^ Xg, 

\A^;^ - Xi\ > kS. 

Pick an "H-orthonormal basis {ve+k}^=o eigenspace, QH, and write 

u (z QTi as u = X^A^o ^kVe+k- Then since is an orthonormal basis for P„7i 

n 

i=i 

29 



Then by (KM 



A J — Aj 7^A^ 

Herein X^a =A{ denotes the sum over those indices j from 1 to n for which Xj 
and similarly for the symbol X^a y^A^- Thus for any normalized u E QTi 



Aj— A^ Aj— A^ 



<(l + ^^ \\u-Pnu\\. (A.3) 

This shows that u € QTi. can be approximated by a linear combination of Rayleigh- 
Ritz eigenvectors. 

To proceed further we first assume that Xi is a simple eigenvalue of A with 
normalized eigenvector u = ug. What then needs to be estimated is the difference 
between {PnUg^uf''^) and 1 as n ^ oo. For this purpose note that we are free to 
orient the unit vectors ui and m^"'' so that {P„Ui, u^"^) is real and nonnegative, which 
minimizes \\ug — (P„Mf , |p. Note that this may make ui depend on n by a 

scalar of magnitude 1. Then, since {PnUi,u^[^^) > 

\\ui\\ - (P„u,,4"V^"^ll < \\{PnUe,u'['^)u'['^\\ 

< |lu,|l + |lu,-(P„w,,4"Vi"^ll (A.4) 

is the same as 

|(P„w,,4"Vl| < ||w£-(P„u,,4"V^"^ll, 

so by (|X3l) 

\\ue - <\\w- (P„u,,4"V(")|| + II ((P„u,,4")) - l) 4"^|| 

<2(^1 + -^^ ||u£-P„M,||. (A.5) 

Now again let = A^+i — . . . — Xi^K have multiplicity K + 1 and let r„ be the 
{K + 1) X {K + 1) matrix given by 



I o (") 
{PnVt,+k,u}^, 



k.m—Q 



30 



Further let V and iJ^"-' be the cohimn vectors 



Vl+l 




r (") 1 

(«) 


ve+K 




(") 



V 



The estimate (|A.3p shows that for each k = 0,1, . . . , K, 

K 



El n (") \ (") 



m=0 



m 7i as n — > oo, i.e. 



V 



in H^'^^ as n — > oo. It foUows that for large n, r„ has rank K + 1, and so is invertible. 
Inversion of r„ would yield a suitable basis for QTL to complete the proof provided 
that it is shown that the inverse is bounded with a bound independent of n. It appears 
that a diagonalization procedure which enables us to mimic the simple eigenvalue case 
yields a more direct route to an explicit estimate. So let u = X^aLo '^kVf.+k & QH 
with \\u\\ = 1. Then 



K n K 

PnU = ^ CkPnVl+k = C.kPnVl+k,U-P^)uf 
j = l k=0 



k=Q 



Since r„ has rank K+1, for each fixed j — 
equations in i^T + 1 unknowns, 



1, ...,£ + K the system of K 



J2 Ck{PnVi+k,u'f^) =0, £<t^j<e + K, 



fe=0 



has a one parameter family of solutions, 



For each j ~ £, £ + 1, . . . , £ + K , fix the solution by requiring that J2k=o I*- 
and that the unit vector 



K iJ") 12 



j,i+k I 



K 



E(n) 
Cj.e+k'"e+k 



fc=0 



in QTi satisfies 



We now replace the arbitrarily chosen basis {vjYjt^^ for Q Ti. with the basis {%,„} 
In place of r„ and V we now have 



£+K 

3=e ■ 



A. 



{PnUe+km,u}' 



e+m' 



K 



k,m—Q 
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and 



Ui+l.n 



Ue+K,n 



respectively. Note that by construction 



In (") \ (") 



{PnUl+K,m U. 




This diagonahzation yields the estimate 




— u 



'i+k 



II <2 1 + 




ll"<'+l,n 



PnUe+k,n\\, k = 0,l,...,K, 



exactly as in the one dimensional case, (jA.4p . (jA.Sp . Hence, for each j ~ 1, . . . ,£ there 
exists a normalized eigenvector Uj,n of A such that 
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